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microstructural  features  and  processes,  at  different  length-scales,  to  the  macroscopic-level  ballistic -penetration  resistance  of  PPTA-based 
fabric  or  PPTA-fiber-reinforced  polymer-matrix  composites.  Specifically,  the  role  of  various  material-synthesis-/fiber-processing-induced 
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material  constitutive  model,  the  role  of  the  material  microstructure  (including  defects)  at  different  length-scales  must  be  taken  into  account. 
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1.  Subject  Materials 

The  work  carried  out  under  the  present  grant  involves  high  specific-strength,  high 
specific-stiffness  p-phenylene  terephthalamide  (PPTA)  polymeric  fiber/filament  (e.g. 
Kevlar®,  Twaron®,  etc.)  based  materials  and  structures.  These  materials/structures  are 
commonly  used  in  various  protective  systems  whose  main  requirement  is  a  high  level  of 
penetration  resistance  against  high-kinetic-energy  projectiles  (e.g.  bullets,  mine,  lED  or 
turbine  fragments,  etc.).  Their  high  mass-efficiency  (i.e.  mass-normalized  performance) 
makes  these  materials  and  structures  particularly  suitable  for  use  in  applications  such  as 
protective  garments  for  personnel  extremity  protection,  interior  spall  liners  in  infantry 
vehicles,  and  a  lining/shroud  for  turbine-fragment  containment.  The  polymeric 
fibers/filaments  are  normally  used  as  either  thread  constituents  in  two-dimensional  or  three- 
dimensional  woven-fabric  structures  or  as  reinforcements  in  high-performance  (typically) 
polymer-matrix  composites. 

2.  Problem  Statement  and  Objeetives 

Development  of  the  aforementioned  protection  systems  is  traditionally  carried  out 
using  legacy  knowledge  and  extensive  fabricate-and-test  procedures.  Since  this  approach  is 
not  only  economically  unattractive,  but  is  often  associated  with  significantly  longer  lead 
times,  it  has  gradually  become  complemented  by  the  appropriate  cost-  and  time-efficient 
Computer  Aided  Engineering  (CAE)  analyses.  This  trend  has  been  accelerated  by  the  recent 
developments  in  the  numerical  modeling  of  transient  non-linear  dynamics  phenomena  such 
as  those  accompanying  blast  and  ballistic  loading  conditions.  However,  the  tools  used  in 
these  analyses  themselves  suffer  from  a  number  of  deficiencies/limitations  which  prevent 
these  analyses  from  being  more  widely  utilized.  In  the  context  of  the  use  of  CAE  analyses  for 
development  of  the  aforementioned  protection  systems,  it  is  well  established  that  one  of  the 
main  deficiencies  stems  from  the  inability  of  currently  available  material  models  to 
realistically  represent  the  response  of  these  materials  under  high-deformation  rate,  large- 
strain,  high-pressure  loading  conditions,  the  conditions  typically  encountered  during 
projectile  impact  events.  One  of  the  reasons  for  the  indicated  short  comings  of  the  present 
material  models  is  a  lack  of  inclusion  of  the  contribution  of  various  phenomena  and 
processes  occurring  at  different  length  scales  to  the  overall  behavior/performance  of  the 
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material.  The  main  objeetive  of  the  present  work  was  to  identify  and  quantify  the 
eontributions  of  the  key  molecule-/fibril-/fiber-/yarn-/ply-seale  mierostruetural  features 
(ineluding  the  defeet  strueture)  to  the  meehanieal  performance  of  the  materials  at  hand. 

3.  Approach 

To  address  the  problem  described  in  the  previous  section,  a  multi-length-scale 
computational  approach  combining:  (a)  quantum-mechanical  calculations;  (b)  all-atom 
molecular-level  techniques;  (c)  coarse-grained  meso-scale  analyses;  (d)  discrete-element 
methods;  and  (e)  continuum- level  approaches  are  utilized.  The  computational  methods  and 
tools  utilized  at  different  length-scales  are  combined  into  a  multi-scale  computational 
framework,  which  enables  interactions  and  linking  of  (i.e.  data-exchange  between)  different 
analyses. 

4.  Relevance  to  the  Army 

In  order  to  respond  to  the  new  enemy  threats  and  warfare  tactics,  military  systems,  in 
particular  those  supporting  the  U.S.  ground  forces,  are  being  continuously  transformed  to 
become  faster,  more  agile,  and  more  mobile  so  that  they  can  be  quickly  transported  to 
operations  conducted  throughout  the  world.  Consequently,  an  increased  emphasis  is  being 
placed  on  the  development  of  improved  lightweight  body-armor  and  lightweight  vehicle- 
armor  systems  as  well  as  on  the  development  of  new  high-performance  armor  materials. 
High-performance  fiber-based  materials,  such  as  PPTA-based  materials,  have  been  exploited 
for  both  body-armor  (e.g.  as  soft,  flexible  fiber  mats  for  personal-armor  vests)  and  for  the 
vehicle-armor  systems  (e.g.  as  reinforcements  in  rigid  polymer  matrix  composites,  PMCs,  for 
lightweight  vehicle-armor  systems).  Establishment  of  the  basic  functional  relationships 
between  the  material  micro  structure  (at  different  length-scales  and  the  material’s  ballistic- 
impact  resistance,  the  subject  of  the  present  work,  is  critical  as  the  Army  is  trying  to  develop 
a  new  generation  of  lightweight  protective  materials. 

5.  Accomplishments 

Significant  progress/accomplishments  have  been  achieved  in  the  following  aspects  of 
the  project:  (a)  identification  of  the  key  mierostruetural  length-scales  in  PPTA-based 
materials;  (b)  identification  of  the  key  synthesis-Zprocessing-induced  defects;  (c)  effect  of 
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synthesis-Zprocessing-induced  defects  on  PPTA-fiber  properties;  (d)  effect  of  fiber-Zyam- 
handling-induced  defects  on  PPTA-fiber  properties;  (e)  development  of  a  continuum-level 
physically-based  material  constitutive  model. 

5.1  Identification  of  the  Key  Microstructural  Length-Scales  in  PPTA-Based  Materials 

In  our  work  reported  in  Refs.  [1,  2],  the  hierarchical  nature  of  the  PPTA-based 
materials  was  investigated.  This  work  identified  eight  microstructural  length-scales, 
schematics  and  explanations  of  which  are  provided  in  Figure  1.  The  first  column  in  this 
figure  shows  a  set  of  simple  schematics  of  the  material  microstructure/architecture  at  a  given 
length-scale  along  with  the  labels  used  to  denote  the  main  microstructural  constituents.  In  the 
second  column,  a  brief  description  is  provided  of  the  material  models  used  to  capture  the 
material  behavior  at  the  length-scale  in  question.  A  brief  description  of  the  material 
micro  structure/  architecture  and  the  corresponding  material  models  at  each  of  the  length- 
scales  listed  in  Figure  1  is  provided  below: 

Laminate  length-scale:  At  this  length-scale,  the  material  possesses  no  discemable 
microstructural  features,  i.e.  it  is  completely  homogenized.  An  example  of  the  laminate- 
length  scale  composite  material  model  can  be  found  in  Ref  [3]. 

Stacked-lamina  length-scale:  At  this  length-scale,  the  presence  of  discrete  stacked  laminae  is 
recognized  while  the  material  within  each  lamina  as  well  as  inter-lamina  boundaries  are  kept 
featureless/homogenized.  An  example  of  the  composite  material  model  at  this  length-scale 
can  be  found  in  Ref  [4]. 

Single-lamina  length-scale:  In  this  case,  the  microstructure/architecture  of  the  reinforcement 
and  matrix  phases  are  explicitly  taken  into  account  while  the  two  associated  materials  are 
treated  as  featureless/homogenized.  The  resulting  '"two-phase'’^  composite-lamina  model  for 
each  lamina  is  then  combined  with  a  fully-homogenized  lamina/lamina  interface  model  to 
form  a  stacked-lamina  composite  laminate  structure.  An  example  of  this  type  of  material 
model  can  be  found  in  Ref  [5]. 

Fabric  Unit-cell  length-scale:  At  this  length-scale,  a  closer  look  is  given  to  the  architecture 
of  the  woven  fabric.  Specifically,  details  of  yarn  weaving  and  crimping,  yam  cross-section 
change,  and  yam  sliding  at  the  warp-yam/weft-yam  cross-over  points  are  taken  into  account. 
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Figure  1.  Various  length-scales  and  the  associated  material  model  assumptions/simplifications  used  in 
the  study  of  polymer-matrix  composite  materials  with  high-performance  fiber-based  structures. 
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Yarn  Length-scale 

Each  yarn  is  considered  as  an  assembly 
of  fully-homogenized  mechanically- 
engaged  discrete  fibers 
Fiber-fiber  contact  and  sliding  are 
modeled  explicitly 


^Fiber-axis 
'-^Direction 


Fiber  Length-scale 

Fibers  are  considered  as  an  assembly  of 
nearly-parallel  fibrils  held  together  by 
weak  van  der  Walls  forces 
Adjacent  fibrils  differ  with  respect  to  the 
orientation  of  the  stacked  sheets. 


Fibril  Length-scale 

Fibrils  are  considered  as  crystalline 
thread-like  entities  containing  parallel 
sheets  of  aligned  and  hydrogen 
bonded  molecular  chains. 

Inter-sheet  bonding  is  of  a  van  der  Walls 
character. 


Molecular-level  Length-scale 

•  Individual  molecules  are  represented  as 

assemblies  of  covalently  bonded 
atoms. 

•  The  effect  of  hydrogen  bonding  on  the 

molecular  topology  (e.g.  the  formation 
of  helical  structure  is  considered 
explicitly). 


Figure  1.  continued. 
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Each  yarn  is  fully  homogenized  in  this  case  as  are  the  matrix  and  the  lamina/lamina 
interfaees.  An  example  of  the  model  of  this  type  ean  be  found  in  Refs.  [6-8]. 

Yarn  length-scale:  In  this  case,  the  internal  structure/arehitecture  of  each  yam  is  accounted 
for  explicitly.  In  other  words,  yarns  are  considered  as  assemblies  of  nearly  parallel 
fibers/filaments  whieh  are  meehanically  engaged  by  either  the  applieation  of  a  light  twist  to 
the  yarn  or  by  wrapping  a  fiber  around  the  fiber/filament  assembly.  Typieally,  the  detailed 
mierostmcture  of  the  yarns  is  not  ineorporated  into  material  models  used  in  large-seale 
simulations  of  projectile/protection-system  impaet  interaetions  (due  to  the  unmanageably 
large  eomputer  resourees  needed).  Instead,  eomputational  results  of  the  meehanieal  response 
of  the  individual  yarns,  when  subjeeted  to  a  variety  of  loading  eonditions,  are  used  to  enrieh 
the  homogenized  material  model  for  the  yams.  The  latter  is  subsequently  used  at  the  fabrie 
unit-eell  length  seale.  An  example  of  this  approaeh  ean  be  found  in  Ref  [9]. 

Fiber  length-scale:  At  this  length-scale,  fibers  are  eonsidered  as  assemblies  of  aligned  long- 
ehain  moleeules  whieh  are  held  together  by  non-bond  (van  der  Waals  or  Coulomb)  forees.  In 
the  ease  of  Kevlar®  type  fibers,  there  is  a  substantial  experimental  support  for  the  existenee 
of  fibrils  within  the  fibers.  Fibrils  are  smaller  bundles  of  moleeular  ehains  within  whieh 
ehain  moleeules  are  tightly  bonded  into  a  perfeet  or  nearly  perfeet  erystalline  phase.  Thus,  in 
the  ease  of  Kevlar®  type  fibers,  fibers  can  be  eonsidered  as  an  assembly  of  fibrils.  At  this 
length- seale,  the  material  is  modeled  using  an  all-atom/moleeular  approaeh  within  which  the 
eonstituent  diserete  partieles  (atoms  or  ions)  internet  via  valenee-bond  and  non-bond  forces. 
As  in  the  previous  ease,  the  knowledge  gained  at  the  fiber  length-scale  is  used  to  enrieh  the 
material  description  at  the  next  level  length-scale  (the  yam  length-seale,  in  the  present  ease) 
and  is  not  directly  used  in  the  large-scale  modeling  of  projeetile/protective-stmeture 
interaetions. 

Fibril-level  length-scale:  As  mentioned  above,  within  the  fibrils  the  moleeular  stmeture  is 
erystalline  or  nearly  erystalline.  However,  the  material  at  this  length-seale  (as  well  as  the 
fiber  length-seale)  may  eontain  a  variety  of  material-synthesis-/fiber-manufaeturing-indueed 
mierostmctural  and  topologieal  defects  and  chemical  impurities  whieh  may  signifieantly  alter 
its  properties.  As  in  the  fiber  length-seale  case,  the  material  is  treated  as  a  eollection  of 
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discrete  interaeting/bonded  partieles  and  analyzed  using  the  aforementioned 
atomie/moleeular  modeling  tools/proeedures. 

Molecular  chain-level  length-scale:  At  this  length-scale,  chemieal  structure  and 
eonformation  of  the  individual  molecules  eonstituting  the  ehain  is  analyzed  using  the 
aforementioned  atomie-Zmolecular- modeling  tools/proeedures.  The  main  goal  of  the  analysis 
at  this  length-scale  is  to  identify  the  most-likely  moleeular  eonformations  present  in  the 
moleeules  and,  thus,  fibrils.  This  greatly  reduees  the  computational  cost  expended  at  the 
fibril  and  fiber  length-scales. 

5.2  Identification  of  the  Key  Synthesis-ZProcessing-Induced  Defects 

As  in  most  engineering  materials,  properties  of  PPTA  fibers  are  greatly  affected  by 
the  presenee  of  various  defeets/flaws.  These  flaws  ean  be  introdueed  during  PPTA  synthesis 
and  PPTA-fiber  fabrieation  proeesses.  In  the  present  work,  the  key  steps  assoeiated  with 
PPTA  synthesis  and  PPTA-fiber  processing  have  been  analyzed  in  detail.  This  analysis 
identified  the  main  synthesis-Zproeessing-indueed  defeets  in  PPTA  fibers.  A  summary  of  the 
PPTA  fiber  most  eommon  defeets,  their  dimensionality,  their  eause,  ways  of  redueing  their 
density,  and  their  typieal  eoneentrations  is  provided  in  Table  1.  These  defeets  are  expeeted  to 
have  a  profound  effeet  on  the  fiber  properties,  as  well  as  on  the  properties  of  eoarser-seale 
materials  and  struetures  eontaining  PPTA  fibers  (e.g.  yams,  fabrie-plies,  lamina  and 
laminates). 

5.3.  Effect  of  Synthesis-ZProcessing-Induced  Defects  on  PPTA-Fiher  Properties 

In  this  section,  a  brief  summary  is  provided  of  the  main  results  obtained  in  the  present 
work,  in  which  all-atom  molecular-level  modeling  was  used  to  establish  the  effect  of  various 
defects  on  the  material  properties  (primarily  strength  and  stiffness).  This  involved  the  use  of 
a  procedure  to  establish  relationships  between  topology/kinematics  and  properties  of  PPTA 
fibrils,  fibers  and  yams.  Figures  l(e)-(g). 

Within  the  analysis  carried  out  in  the  present  work,  it  was  demonstrated  that  the 
longitudinal  stiffness,  strength  and  ductility  of  the  fibrils  make  the  dominant  contribution  to 
the  corresponding  fiber  and,  in  turn,  yarn  properties.  As  far  as  the  transverse  strength  of 
fibers  and  yarns  are  concerned,  they  are  found  to  be  controlled  by  weak  inter- fibril  (and  in 
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Table  1.  A  summary  of  most  common  defects  found  in  Kevlar®-based  materials 


Flaw  Group 

Flaw 

Causes 

How  causes  can  be 
remediedfmitigated 

Range  for  concentration 

■COOH 

H2SO4  catalyzed 
hydrolysis  causing  PPTA 
chain  scission.  Na+ 
deficiency  with  respect  to 
complete  neutralization  of 
sidefend  acidic  groups. 

Use  concentrated  H2S04for 
dope  preparation.  Shorten  the 
fiber  wash  time 

0.35  per  PPTA  chain  for 
each  flaw* 

(-350  ppm-mass-based) 

Isolated  chain 
ends 

(point  defect) 

-NH2 

H2SO4  catalyzed 
hydrolysis  causing  PPTA 
chain  scission.  Na+ 
deficiency  with  respect  to 
complete  neutralization  of 
side/end  acidic  groups. 

Use  higher  concentration  NaOH 
solution 

0.35  per  PPTA  chain  for 
each  flaw* 

(-350  ppm-mass-based) 

■COONa" 

COOH  neutralization  with 
Na+ 

No  remedy  required  since  this  is 
one  of  the  preferred  chain  ends 

1.1  per  PPTA  chain* 
(-1100  ppm-mass-based) 

-NH3^HS04- 

Sulfonation  of  the  NH2 
chain  ends 

Increase  the  H2SO4  removal  and 
neutralization  rate 

0.2  per  PPTA  chain* 
(-200  ppm-mass-based) 

Side  Groups 
(point  defect) 

■SO3H 

Exposure  of  PPTA  in  the 
dope  to  concentrated 
H2SO4  (sulfonation) 

Reduce  the  H2SO4  concentration 
in  the  dope 

-1300  ppm  (mass-based) 

■SOrNa* 

Neutralization  of  sulfonic 
acid  side  groups  by  NaOH 

Remedy  may  not  be  required 
since  this  side  group  improves 
fiber  longevity.  However 
mechanical  performance  may  be 
compromised 

-2500  ppm  (mass-based) 

Voids  and 
Interstitials 

Microvoids 

Swelling  induced  by 
hydration  of  intra-fibrillar 
Na2S04 

Increase  the  extent  of  sodium 
salt  dissolution  by  prolonged 
exposure  of  fibers  to  boiling 
water 

-150  ppm  (mass-based) 

(point  defects) 

Mobile  Trapped  H2SO4 

Non-neutralized  or 
unwashed  intra-fibrillar 
H2SO4 

Thorough  washing  in  hot 
solvent  aqueous  bath 

-70  ppm  (mass-based) 

Defect  bands 
(planar  defects) 

NH3"HS04- 
agglomerated  chain 
ends 

Coulombic  forces  induced 
clustering  of  ion- 
terminated  chain  ends 

The  phenomenon  is  not  well 
understood  so  no  remedy  is 
obvious 

One  band  every  40-60nm 
of  fibril  (ca.  3000ppm- 
mass-based)) 

*  Extruded  fibers 
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the  case  of  yarns,  weak  inter- 
fiber)  van  der  Waals  / 
mechanical  interlocking 
forces  and  to  be  less 
dependent  on  the  fibril 
transverse  mechanical 

properties. 

The  results  obtained 
in  this  portion  of  the  work 
can  be  summarized  as 
follows: 

(a)  fibril  stiffness  was 
found  to  be  affected  by  both 
the  type  and  concentration  of 
the  defects.  An  example  of 
the  typical  results  pertaining 
to  the  effect  of  defect 
type/concentration  on  the 
PPTA-fibril  stiffness  is 
depicted  in  Figures  2(a)-(b). 
In  Figure  2(a),  the  case  of- 
COOH  chain  end  defect  was 
considered,  while  in  Figure 
2(b)  the  case  of  -SO3H  side 
group  defect  was  analyzed; 

(b)  As  far  as  the  fibril 
axial  strength  was 
concerned,  it  was  found  to 

be  affected  by  defect  type  and  the  size  and  composition  of  the  largest  defect  cluster. 
Examples  of  typical  results  obtained  are  shown  in  Tables  2  and  3.  The  clusters  composed  of 
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Table  2.  A  summary  of  the  effect  of  the  most  common  defects  found  in  Kevlar®-based  materials  and  their 
(same  defect-type)  clusters  on  the  PPTA  fibril  axial  tensile  strength  (in  GPa) 


Flaw  Group 

Flaw 

Number  of  defects  in  the  largest  cluster 

1 

2 

3 

Isolated  chain  ends 
(point  defect) 

-COOH 

6.0 

5.6 

5.0 

-COONa+ 

-NH2 

6.1 

5.8 

5.2 

-NH3+HS04- 

Side  Groups 
(point  defect) 

-SO3H 

5.8 

5.1 

4.1 

-S03Na+ 

Voids  and 
Interstitials 
(point  defects) 

Microvoids 

3.8 

N/A 

N/A 

Mobile  Trapped 
H2SO4 

3.9 

N/A 

N/A 

Defect  bands 
(planar  defects) 

NH3"HS04- 
agglomerated  chain 
ends 

3.2 

N/A 

N/A 

Table  3.  A  summary  of  the  effect  of  the  composition  and  the  size  of  the  mixed  S03H/S03'Nafi3ased  side 
group  and  COOH/COO'Na^  based  chain  end  defect  clusters  on  the  PPTA  fibril  axial  tensile  strength  (in 
GPa) 


Number  of  SOsH/SOs  Na+defects 
in  the  cluster 

1 

2 

3 

Number  of 
COOH/COONa+ 
defects  in  the 
cluster 

1 

5.5 

4.8 

3.8 

2 

5.3 

4.4 

3.4 

3 

4.5 

3.8 

2.7 
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defects  of  the  same  type  are 
considered  in  Table  2  while  the 
effect  of  mixed-type  defect 
(chain-end  and  side-group) 
clusters  is  considered  in  Table 
3.  Examination  of  the  results 
reported  in  Tables  2  and  3 
reveals  that  the  extent  of  fibril 
axial-strength  loss  increases 
with  the  cluster  size/strength; 

(c)  The  aforementioned 
finding  that  the  fibril  axial 
strength  is  a  function  of  the 
largest  defect  cluster  and  its 
composition  indicated  that 
fibril  strength  is  a  stochastic 
property.  Elence,  in  the  present 
work  results  like  the  ones 
displayed  in  Tables  2  and  3  are 
combined  with  the  prototypical 
defect  concentration  densities, 
Table  1,  in  order  to  derive  the 
associated  PPTA  fibril- strength 
distribution  function,  Figures 
3(a)-(b); 

(d)  In  accordance  with 
the  established  topological 

relationships  between  the  fibrils,  fibers  and  yams,  PPTA  fibers  are  found  to  be  transversely 
isotropic  (with  the  fiber  axis  being  the  unique  material  direction)  although  the  constituent 
fibrils  possess  orthotropic  symmetry.  Consequently,  (i)  fiber  axial  stiffness  and  strength, 
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and  (Ti  are  found  to  be  nearly  equal  to  their  fibril  eounterparts;  (ii)  the  fiber  transverse  normal 
stiffness,  £22=^^^,  is  found  to  be  nearly  equal  to  the  more  eompliant  £’33  fibril  modulus;  and 
(iii)  All  three  shear  moduli  G^,  and  G23  are  found  to  be  nearly  equal  to  the  same  value, 
the  value  whieh  is  eontrolled  by  the  eompliant  PPTA  inter-sheet  sliding  resistanee  (i.e.  G12  of 
the  fibrils);  and 

(e)  To  assess  the  yarn-level  stiffness  and  strength  properties  from  the  known  fiber 
properties,  additional  effeets  arising  from  the  fiber-twist-indueed  inter- fiber  frietion  had  to  be 
taken  into  aeeount.  This  proeedure  mainly  affeeted  the  yarn  axial  strength  (or  more  preeisely, 
yarn  axial-strength  distribution).  It  was  generally  found  that  the  fibril-level  axial  strength 
distribution  funetion  displayed  in  Figures  3(a)-(b)  was  modified  in  the  yam  ease  by  an  inter- 
fiber- frietion  dependent  (positive)  strength  shift.  For  example,  for  a  prototypieal  extent  of 
fiber  twist,  an  inerease  in  the  strength  of  ca.  0.2  GPa  was  found. 

5.4  Effect  of  Fiber-/Yarn-Handling-Induced  Defects  on  PPTA-Fiber  Properties 

During  the  weaving  proeess  used  to  produee  the  PPTA  fabrie,  PPTA  fibers/yams  (in 
particular,  the  warp  fibers/yarns)  can  sustain  substantial  damage  due  to  high  bending  and 
contact  stresses.  The  resulting  defects  can  also  degrade  mechanical  properties  of  the  PPTA- 
based  materials,  and  their  effect  has  to  be  accounted  for  in  the  analysis  of  the  ballistic- 
protection  resistance  of  the  PPTA-based  materials.  To  address  this  problem,  a  series  of  all¬ 
atom  and  coarse-grained  molecular-level  computational  analyses  is  carried  out  in  the  present 
work  [11-13].  Specifically,  the  effect  of  the  preludial  handling-induced  deformation  modes 
such  as  fiber/yarn  transverse  compression  due  to  inter-yarn  contacts,  fiber/yam  twist/axial- 
compression,  and  fiber/yam  bending  (encountered  primarily  in  warp  yarns)  and  the 
associated  formation  of  kink-bands  has  been  investigated.  The  results  obtained  revealed: 

(a)  high  contact  stresses  can  induce  topological/crystallographic  defects  and,  in  turn, 
severely  degrade  the  longitudinal-tensile  strength  of  PPTA  fibers/yams,  and  significantly 
alter  the  associated  PPTA- fiber  longitudinal- strength  probability  distribution  function; 

(b)  application  of  the  prior  axial  torsion  to  the  PPTA  imperfect  single-crystalline 
fibrils  degrades  their  longitudinal-tensile  strength  and  only  slightly  modifies  the  associated 
probability  density  function.  This  effect  of  the  preludial  axial  torsion  has  been  linked  to  the 
permanent  changes  in  the  fibril-/fiber-microstmcture  caused  by  this  deformation  mode.  An 
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example  of  some  changes  is 
depicted  in  Figures  4(a)-(b). 
The  results  displayed  in 
Figure  4(a)  pertain  to  the 
microstructure  before 

deformation,  while  Figure 
4(b)  shows  the 

micro  structure  after  the 
application  and  subsequent 
removal  of  the  axial  torsion 
associated  with  a  maximum 
shear  strain  of  0.1. 

Examination  of  the  results 
displayed  in  Figures  4(a)-(b) 
shows  that  the  application  of 
simple  shear  gives  rise  to 
some  permanent  changes  in 
the  PPTA  microstructure. 
Specifically,  under  axial 
torsion,  an  increase  in  the 

extent  of  disorder  is 

observed,  particularly  in  the 
region  adjacent  to  the  single 
defect  present  in  the 

computational  cell.  This  disorder  is  accompanied  by  a  reduction  in  the  number  of  hydrogen 
bonds.  These  residual  changes  in  the  PPTA-crystal  microstructure  are  a  probable  reason  that 
the  prior  axial  torsion  causes  the  changes  in  the  PPTA-fiber  longitudinal-tensile  strength 
probability  distribution  function; 

(c)  The  preludial  axial  compression  can  introduce  (residual)  kink  bands  into  the  PPTA 
fibers  and,  in  turn,  significantly  degrade  fiber  longitudinal-tensile  strength  and  its  statistical 


Figure  4.  Molecular-level  microstructure  in  imperfect  PPTA 
crystal:  (a)  before;  and  (b)  after  the  application  and  removal  of 
axial  torsion  strain  of  0. 1 . 
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distribution.  An  example  of  the  axial-eompression-induced  shear  bands  is  depicted  in  Figure 
5; 

(d)  fiber-fibrillation  can  also  take  place  during  axial  compression,  and  give  rise  to  the 
longitudinal-tensile  strength  and  degradation,  and  changes  in  the  strength  probability  density 
function.  An  example  of  the  axial-compression- induced  fiber- fibrillation  is  shown  in  Figure 
6. 

5.5  Development  of  a  Continuum-Level  Physieally-Based  Material  Constitutive  Model 

PPTA-fiber-reinforced  polymer  matrix  composite  materials  display  quite  complex 
deformation  and  failure  behavior  under  ballistic/blast  impact  loading  conditions.  This 
complexity  is  generally  attributed  to  a  number  of  factors  such  as:  (a)  hierarchical/multi¬ 
length  scale  architecture  of  the  material  micro  structure;  (b)  non-linear,  rate-dependent  and 
often  pressure-sensitive  mechanical  response;  and  (c)  the  interplay  of  various  intrinsic 
phenomena  and  processes  such  as  fiber  twisting,  inter-fiber  friction/sliding,  etc.  Material 
models  generally  employed  in  the  computational  engineering  analyses  of  ballistic-Zblast- 
impact  protective  structures  made  of  this  type  of  material  do  not  typically  include  many  of 
the  aforementioned  aspects  of  the  material  dynamic  behavior.  Consequently,  discrepancies 
are  often  observed  between  computational  predictions  and  their  experimental  counterparts. 
To  address  this  problem,  the  results  of  an  extensive  set  of  molecular-level  computational 
analyses  carried  out  in  the  present  work  and  overviewed  earlier,  regarding  the  role  of  various 
microstructural/morphological  defects  on  the  PPTA  fibril/fiber/yarn  mechanical  properties, 
are  used  to  upgrade  one  of  the  continuum-level  material  models  for  fiber-reinforced 
composites  developed  by  Yen  [14]  of  the  Army  Research  Laboratory,  Aberdeen  Proving 
Ground,  MD.  The  results  obtained  in  the  upgraded  PPTA  continuum-level  material  model 
[14-18]  show  that  the  response  of  the  material  is  significantly  affected  as  a  result  of  the 
incorporation  of  micro  structural  effects  both  under  quasi-static  simple  mechanical  testing 
condition  and  under  dynamic  ballistic-impact  conditions.  This  is  demonstrated  in  Figure  7,  in 
which  typical  results  pertaining  to  the  effect  of  projectile  initial  velocity  on  the  projectile 
residual  velocity  in  the  case  of  the  Yen  [14]  composite  material  model  and  its  present 
modification  are  displayed.  It  is  seen  that  Vso  (defined  here  as  the  projectile  initial  velocity  at 
which  the  projectile  residual  velocity  is  zero)  predicted  by  the  analysis  in  which  the 
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Figure  5.  Coarse-grained  mierostrueture,  during  axial  eompression,  within  the  eomputational  eell 
eontaining  (initially)  defeet-free  virgin  PPTA  fiber  at  the  longitudinal  streteh  values  of:  (a)  0.975;  (b) 
0.950;  (e)  0.925;  and  (d)  0.900. 
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Figure  5.  continued. 
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Figure  6.  Coarse-grained  microstructure,  during  axial  compression,  within  the  computational  cell 
containing  (initially)  defective  virgin  PPTA  fiber  at  the  longitudinal  stretch  values  of:  (a)  0.975;  (b) 
0.950;  (c)  0.925;  and  (d)  0.900.  The  fiber  within  the  computational  cell  contains  a  cluster  of  side- 
group  defects. 
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modified  composite 
material  model  is  used 
is  ca.  5%  lower  than  its 
counterpart  associated 
with  the  use  of  the 
original  Yen  [14] 
composite  material 
model.  This  effect  is 
deemed  as  being 
significant  and  justifies 
the  present  efforts 
aimed  at  incorporating 
fine-scale 
microstructural  effects 
into  the  continuum- 
level  material  model. 

An  example  of 

the  prototypical  results  obtained  in  the  present  work,  pertaining  to  the  ballistic-impact  of  a 
right-circular-cylindrical  rigid  projectile  on  a  PPTA-fiber-reinforced  composite  laminate 
containing  a  poly-vinyl-ester-epoxy  matrix,  is  given  in  Figures  8(a)-(d).  The  results 
displayed  in  Figures  8(a)-(d)  pertain  to  the  spatial  distribution  of  the  composite-laminate 
material  during  impact.  It  should  be  noted  that  in  order  to  enable  monitoring  of  the 
penetration  process,  the  target  plate  is  made  transparent.  The  results  displayed  in  Figures 
8(a)-(d)  clearly  reveal  the  development  of  the  back  face  bulge.  In  general,  it  was  found  that 
details  regarding  the  spatial  distribution  and  the  temporal  evolution  of  the  target  material 
during  impact  are  significantly  affected  by  the  modifications  made  to  the  original  Yen  [14] 
composite  material  model. 
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Figure  8.  Typical  results  pertaining  to  the  temporal  evolution  of  the  spatial  distribution  of  the 
composite-laminate  material  during  impact  by  a  right-circular  cylindrical  rigid  projectile. 


21 


6.  Bibliography 

1.  M.  Grujicic,  W.  C.  Bell,  P.  S.  Glomski,  B.  Pandurangan,  C-F.Yen  and  B.  A. 
Cheeseman,  "^Filament-level  Modeling  of  Aramid-based  High-Performance  Structural 
Materials,”  Journal  of  Materials  Engineering  and  Performanee,  20,  1401-1413,  2011. 

2.  M.  Grujieie,  W.  C.  Bell,  P.  S.  Glomski,  B.  Pandurangan,  C-F.  Yen  and  B.  A. 
Cheeseman,  ""Multi-length  Scale  Computational  Derivation  of  Kevlar®  Yarn-level 
Material  Model,”  Journal  of  Materials  Seience,  46,  4787-4802,  2011. 

3.  M.  Grujieie,  B.  Pandurangan,  K.  L.  Koudela  and  B.  A.  Cheeseman,  ""A 
Computational  Analysis  of  the  Ballistic  Performance  of  Light-Weight  Hybrid- 
Composite  Armor,”  Applied  Surfaee  Seienee,  253,  730-745,  2006. 

4.  M.  Grujieie,  T.  Fie,  H.  Marvi,  B.  A.  Cheeseman,  C.-F.  Yen,  ""A  Comparative 
Investigation  of  the  Use  of  Laminate-level  Meso-scale  and  Fracture-mechanics 
Enriched  Meso-scale  Composite-material  Models  in  Ballistic  Resistance  Analyses,” 
Journal  of  Materials  Seienee,  45,  3136-3150,  2010. 

5.  M.  Grujieie,  D.  C.  Angstadt,  Y.-P.  Sun  and  K.  L.  Koudela,  ""Micro-mechanics 
Based  Derivation  of  the  Materials  Constitutive  Relations  for  Carbon  Nanotube 
Reinforced  Poly-Vinyl-Ester-Epoxy  Based  Composites,”  Journal  of  Materials  Seience, 
42,  2007,  4609-4623. 

6.  M.  Grujicic,  G.  Arakere,  T.  Fie,  M.  Gogulapati  and  B.  A.  Cheeseman,  ""A 
Numerical  Investigation  of  the  Influence  of  Yarn-Level  Finite-Element  Model  on 
Energy  Absorption  by  a  Flexible-fabric  Armor  During  Ballistic  Impact,”  Journal  of 
Materials:  Design  and  Applications,  222,  2008,  259-276. 

7.  M.  Grujicic,  W.  C.  Bell,  G.  Arakere,  T.  Fie,  X.  Xie  and  B.  A.  Cheeseman, 
""Development  of  a  Meso-scale  Material  Model  for  Ballistic  Fabric  and  its  Use  in 
Flexible- armor  Protection  Systems,”  Journal  of  Materials  Engineering  and 
Performance,  19-1,  2010,  22-39. 

8.  M.  Grujicic,  W.  C.  Bell,  T.  He  and  B.  A.  Cheeseman,  ""Development  and 
Verification  of  a  Meso-scale  based  Dynamic  Material  Model  for  Plain-woven  Single- 
ply  Ballistic  Fabric,”  Journal  of  Material  Science,  43,  2008,  6301-6323. 

9.  B.  A.  Cheeseman,  C.  F.  Yen,  B.R.  Scott,  B.  Powers,  T.A.  Bogetti,  B.  LaMatina,  Y. 
Duan,  M.  Keefe,  Y.  Miao  and  Y.  Wang,  ""From  Filaments  to  Fabric  Packs  - 
Simulating  the  Performance  of  Textile  Protection  Systems,”  ARE  Report 
ADM002075,  2006. 

10.  S.  M.  Lee,  Handbook  of  Composite  Reinforcements,  John  Wiley  and  Sons,  ISBN: 
978-0-471-18861-2,  1993. 

11.  M.  Grujicic,  S.  Ramaswami,  J.  S.  Snipes,  R.  Yavari,  G.  C.  Lickfield,  C.-F.  Yen 
and  B.  A.  Cheeseman,  ""Molecular-Level  Computational  Investigation  of  Mechanical 
Transverse  Behavior  of  p-phenylene  terephthalamide  (PPTA)  Fibers,”  Multidiscipline 
Modeling  in  Materials  and  Structures,  9,  462-498,  2013.  DOI  10.1108/MMMS-ll- 
2012-0018 


22 


12.  M.  Grujicic,  R.  Yavari,  J.  S.  Snipes,  S.  Ramaswami,  C.-F.Yen,  and  B.  A. 
Cheeseman,  "^Molecular-Level  Study  of  the  Effect  of  Prior  Axial  Compression/Torsion 
on  the  Axial-Tensile  Strength  ofPPTA  Fibers  f  Journal  of  Materials  Engineering  and 
Performance,  22,  3269-3287,  2013.  DOI:  10. 1007/sl  1665-013-0648-2 

13.  M.  Grujicic,  S.  Ramaswami,  J.  S.  Snipes,  R.  Yavari,  C.-F.  Yen  and  B.  A. 
Cheeseman,  ""Axial-Compressive  Behavior,  Including  Kink-Band  Formation  and 
Propagation,  of  Single  p-Phenylene  Terephthalamide  (PPTA)  Fibers  f  Advances  in 
Materials  Science  and  Engineering,  Volume  2013,  Article  ID  329549,  2013.  DOE 
10.1155/2013/329549 

14.  C.  F.  Yen,  ""A  Ballistic  Material  Model  for  Continuous-Fiber  Reinforced 
Composites  f  International  Journal  of  Impact  Engineering,  46,  1 1-22,  2012. 

15.  M.  Grujicic,  B.  Pandurangan,  J.  S.  Snipes,  C-F.Yen  and  B.  A.  Cheeseman,  ""Multi- 
Length  Scale  Enriched  Continuum-Level  Material  Model  for  Kevlar® -Fiber 
Reinforced  Polymer-Matrix  Composites  f  Journal  of  Materials  Engineering  and 
Performance,  22,  681-695,  2013.  DOE  10.1007/sl  1665-012-0329-6 

16.  M.  Grujicic,  A.  Hariharan,  B.  Pandurangan,  C-.F.  Yen,  B.  A.  Cheeseman,  Y. 
Wang,  Y.  Miao  and  J.  Q.  Zheng,  ""Fiber-level  Modeling  of  Dynamic  Strength  of 
Kevlar®  KM2  Ballistic  Fabric  f  Journal  of  Materials  Engineering  and  Performance, 
21,  1107-1119,  2012.  DOE  10.1007/sl  1665-01 1-0006-1. 

17.  M.  Grujicic,  B.  Pandurangan,  W.  C.  Bell,  C-F.Yen  and  B.  A.  Cheeseman, 
""Application  of  A  Dynamic-mixture  Shock-wave  Model  to  the  Metal-matrix 
Composite  Materials  f  Materials  Science  and  Engineering  A,  528,  8187-8197,  2011. 

18.  M.  Grujicic,  W.  C.  Bell,  B.  Pandurangan,  C-.F.  Yen  and  B.  A.  Cheeseman, 
""Computational  Investigation  of  Structured  Shocks  in  Al/SiC-particulates  Metal 
Matrix  Composites  f  Multidiscipline  Modeling  in  Materials  and  Structures,  7,  469- 
497,  2011. 


23 


